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ABSTRACT 



Aims. To investigate the star formation activity and the gas and stellar dynamics on scales of a few parsecs in the nucleus of 
the Circinus Galaxy. 

Methods. Using the adaptive optics near infrared integral field spectrometer SINFONI on the VLT, we have obtained 
observations of the Circinus galaxy on scales of a few parsecs and at a spectral resolution of 70kms~^ FWHM. The physical 
\^^ ' properties of the nucleus are analyzed by means of line and velocity maps extracted from the SINFONI datacube. Starburst 
f^ models are constrained using the Br7 flux, stellar continuum (as traced via the CO absorption bandheads longward of 2.3/im), 
and radio continuum. 

Results. The similarity of the morphologies of the H2 1-0 S(l) 2.12^m and Br7 2.17/im lines to the stellar continuum 
and also their kinematics, suggest a common origin in star formation. Within 8pc of the AGN we find there has been 
a recent starburst in the last 100 Myr, which currently accounts for 1.4% of the galaxy's bolometric luminosity. The 
similarity of the spatial scales over which the stars and gas exist indicates that this star formation is occuring within 
the torus; and comparison of the gas column density through the torus to the maximum possible optical depth to the 
stars implies the torus is a clumpy medium. The coronal lines show asymmetric profiles with a spatially compact narrow 
component and a spatially extended blue wing. These characteristics are consistent with some of the emission arising in 
clouds gravitationally bound to the AGN, and some outflowing in cloudlets which have been eroded away from the bound clouds. 
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1. Introduction are exactly those that can be resolved with SINFONI in 

the nearest AGN. 
In the context of active galactic nuclei, star formation ac- 
tivity and the gas and stellar dynamics on scales of a few The Circinus galax y, at a distance of 4.2 ± 0.8 Mpc 
parsecs to a few tens of parsecs can be counted among the (Freeman et al. | 1977| ), is an ideal subject to study be- 
main debated issues. The unified model of active galaxies cause of its proximity (1" = 20 pc) I t is a large, highly 
(see Lawrence USHZI for a review) assumes that the inner inclined {i = 65° Freeman et al. 119771), spiral galaxy that 
region of Seyfert 2 galaxies is comprised of a dense circum- ^^osts both a typical Seyfert 2 nucleus and a circum nuclea r 



nuclear torus that hides the nucleus and the broad-line starburst on scales of 100-200 pc (Maiolino et al. 119981) 

region from our line of sight at near-infrared and optical Evidence for an obscured Seyfert 1 nucleus is provided 

wavelengths. The size scales on which models predict the by the finding of a broad (FWHM = 330 kms i) Ha line 

canonical torus vary from an inner edge at 1 pc out to component in polarized light (Oliva et al.^MB- Their pic- 

sever al tens of parsecs (P ier fc Krolik .1992. Nenkova et al. tnre is also supported by recent X-ray observations. Those 

\Mm Schartmann et al. ^UU^. These crucial size scaks above 10 keV suggest direct X-ray detection of the nucleus 

through a column density of 4 x 10^'' cm ^ (Matt et al. 

Send offprint requests to: F. Mueller Sanchez UJMM- The X-ray spectrum below 10 keV exhibits a flat 

e-maU: frankmueller@mpe.mpg.de continuum and a very prominent iron fine, indicative of 

* Based on observations collected at the European Southern Compton scattering and fluorescent emission from gas il- 

Observatory, Chile luminated by an obscured X-ray continuum source (Matt 
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et al. I1996|l . Circinus shows highly ionized gas extending 
along the minor axis of the galaxy, with a morphology that 
is reminiscent of the ionization cones seen in other Seyfert 
galaxies (Marconi et al. 2994a). H2O maser emission has 
been detected, and the masing gas traces a thin accretion 
disk about 0.4 pc in radius, with, in addition, a fraction of 
the masers originating outside the disk, in what appears 
to be an outflow within ~1 pc of the nucleus and aligned 
with the ionization cone (Greenhill et al. l2(JU3fl . Optical 
and near-infrared spectrophotometry of the nucleus show 
a typical Seyfert spectrum, including strong coronal lines 
(Oliva et al. lHMl Prieto et al. lMI^ . There are also lines 
from H2 and low-excitation ionic species, both believed 
to be associated with starforming regions. The distribu- 
tion and kinematics of the Br7 line have been interpreted 
in terms of ongoing star formation activity within a few 
tens of parsecs of the active nucleus (Maiohno et al.|T998). 
A young stellar population with an age between 4 x 10^ 
and 1.5 x 10® was found between these scales. Recent ob- 
servations of the Circinus Galaxy in the range between 
1 — 10/J,m (Prieto et al. I2004|l . resolve a Ks-hand source 
with a FWHM of ~ 2 pc and a spectral energy distribution 
compatible with a dust temperature of 300 K. 

This paper presents high-resolution, nearly diffraction- 
limited integral field spectroscopic data of the Circinus 
galaxy in the i^-band observed with the adaptive 
optics (AO) assisted imaging spectrograph SINFONI 
HBonnet et al. 20041 [Eiscnhaucr ct al. 2003b|)- The phys- 
ical scale associated with our spatial resolution is 4pc, 
allowing us to investigate the properties of the nucleus 
and its interaction with the circumnuclear environment, 
highlighting the relationship between star formation and 
galactic nuclear activity. 

We describe our observations and the data reduction 
procedure in Sect. El Sect. 01 describes the features ob- 
served in the spectrum of the galaxy and the methods 
we used to estimate the spatial resolution of our observa- 
tions. In Sect. ^ we analyze the morphology of the nuclear 
source and the origin of the observed K-hand continuum. 
Sect. discusses the morphologies of the emission and the 
absorption, and the age of the starburst at the nucleus. 
Then, in Sect. we discuss the kinematics of the emission 
and the mass distribution in the galaxy. Following these 
results, Sect. deals with the gas density at the nucleus 
and the physical characteristics of the torus. A study of 
the detected coronal hues is presented in Sect.JSl and we 
present our conclusions in Sect. [51 

2. Observations and data reduction 

The data presented here were obtained on 15 Jul 2004 
during commissioning of SINFONI IjBonnet et al. 200'4l 
lEisenhauer et al. 2003b|l on the VLT UT4. The instru- 
ment consists of a cryogenic near infrared integral field 
spectrometer SPIFFI l|Eisenhauer et al. 2003ajl coupled 
to a visible curvature adaptive optics (AO) system 
IjBonnet et al. 2003)l . SINFONI performs imaging spec- 
troscopy by cutting the two-dimensional field of view into 



32 slices and then rearranging each of the slices (slitlets) 
onto a one-dimensional pseudo longslit, which is dispersed 
by a grating wheel and at last the spectra are imaged on 
the detector. 

The AO module was able to correct on the nucleus of 
Circinus (for which it measured R=14.1mag) in seeing of 
^^0.5", to reach a resolution of 0.2" in the iiT-band (see 
Sectional). This performance is good considering that the 
optical nucleus of Circinus is rather extended rather than 
point-like. With the appropriate pixel scale selected, the 
spectrograph was able, in a single shot, to obtain spectra 
covering the whole of the i^-band (approximately 1.95- 
2.45 fj,m) at a spectral resolution of i? ^ 4200 for each 
0.0125" X 0.025" pixel in a 0.80" x 0.80" field of view. A 
total of 6 sky and 6 on-source exposures of 300 sec each 
were combined to make the final data cube with a total 
integration time of 1800 sec. 

The data were reduced using the SINFONI custom re- 
duction package SPEED (Abuter et al. I2005|l . This per- 
forms all the usual steps needed to reduce near infrared 
spectra, but with the additional routines for reconstruct- 
ing the data cube. Following subtraction of the sky frames 
from the on-source frames, the data were flatfielded and 
corrected for dead/hot pixels. The data were then interpo- 
lated to linear wavelength and spatial scales, after which 
the slitlets were aligned and stacked up to create a cube. 
Finally the atmospheric absorption was compensated us- 
ing the AOV star HD 190285. Flux calibration was per- 
formed by comparison to the high spatial resolution broad 
band Kg image and photometry obtained with NAOS- 
CONICA of 11.4 mag in a 0.38"aperture (Prieto et al. 
2004), which was also cross-checked with 2MASS photom- 
etry. 

No additional point-spread function (PSF) calibration 
frames using stars were taken. This is primarily because, 
although in principle one can match the brightness of a 
calibration star on the wavefront sensor to the AGN, it is 
not possible to replicate either the spatial extent of the 
AGN or the background galaxy light associated with it - 
resulting in a potentially considerable mismatch between 
the science and calibration PSFs IjDavies et al. 2004)l . The 
spatial resolution has instead been measured using the 
methods described in Section IXTl 

3. Spectroscopic features and spatial resolution 

The purpose of this section is to present general results 
extracted from the data to provide an overview of the 
features observed in the nuclear region of Circinus. A more 
detailed analysis is given in subsequent sections. 

3.1. Spatial Resolution 

A good way to estimate the spatial resolution is by using a 
broad line in the spectrum, such as a potential observation 
of Br7 emission from the BLR, which at the distance of 
Circinus is expected to be spatially unresolved. This has 
the advantage of measuring the spatial resolution directly 
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Flux*" 




(/im) 


(10-^« Wm-2) 


Ha 1-0 S(3) 


1.9576 


6.0 


[Sivi] 


1.9634 


53.9 


Ha 1-0 S(2) 


2.0338 


3.8 


[Alix] 


2.040 


10.0 


[Hei] 


2.0587 


5.0 


Ha 1-0 S(l) 


2.1218 


9.0 


Br7 


2.1661 


14.1 


Ha 1-0 S(0) 


2.2233 


3.4 


[Caviii] 


2.3213 


36.9 


Ha 1-0 Q(l) 


2.4066 


3.4 



° Wavelengths are in the rest frame 

*" Uncertainities are approximately 3.2 x 10"^'^ Wm~^ 

Table 1. Measured nuclear emission line fluxes for 
Circinus. All measurements are given for a 0.8" circular 
aperture centered on the continuum peak. 



from the science frames, and includes all effects associated 
with shifting and combining the cube. Nevertheless, as can 
be seen in Figure ^ we have detected only narrow lines, 
which are spatially resolved in our observations. As a re- 
sult, and due to the lack of other point sources in the field, 
we used indirect means to derive our spatial resolution, 
making use of a NACO image of the galaxy from Prieto 
et al. H2(JU4|I which contains unresolved star clusters over a 
field of view of 27" x 27". Two approaches were used to es- 
timate the PSF of our observations by means of the NACO 
image. One method of estimating the SINFONI PSF ef- 
fectively consists of deconvolving the SINFONI continuum 
image with the intrinsic image of the galaxy, since: 

IMs = PSFs ® IMintr 

where IMs is the Circinus image from SINFONI, PSFs is 
the point spread function, and IMintr is the intrinsic image 
of the galaxy. We obtained the intrinsic image in a simi- 
lar way, by first deconvolving the NACO image with the 
NACO PSF. In the NACO image, the smooth underly- 
ing galaxy profile was approximated by isophotal analysis 
using the isophote package from IRAF. This was then sub- 
tracted, leaving only the compact sources in the field. The 
narrowest point in this subtracted image was used as the 
NACO PSF, which was measured directly to have 0.14" 
FWHM. The deconvolution of the NACO image with its 
PSF was performed using the Lucy algorithm fLucv [T^71)l 
also implemented in IRAF. This deconvolved image was 
rebinned to our pixel scale of 0.0125"/pix. For the final 
stage, instead of deconvolving IMs with IMintr, we fit the 
SINFONI image with the intrinsic image convolved with 
the parameterized PSFs. This avoids the noise amplifi- 
cation which is an inherent feature of deconvolution. The 
PSFs could be well matched by a symmetrical moffat func- 
tion. Mismatches were minimized by varying the parame- 
ters (center, scale, a and (3) of a moffat function defined 
as I(r) = (l + ( — )^) and then performing a x^ mini- 



mization test, resulting in a FWHM spatial resolution of 
0.20" (- 3.9 pc) with a = 0.148" and [3 = 1.8 for this 
method. An elliptical moffat function was also tested as a 
possible PSF for our observations. However, this did not 
improve significantly the fit and so was discarded. 

In order to avoid the first deconvolution of the method 
described above, the second approach to derive the 
SINFONI PSF consists of convolving the NACO PSF with 
a parameterized degradation function, the latter being ob- 
tained by fitting the SINFONI image with the NACO 
image convolved with the same function. A symmetrical 
moffat function was used as degradation function. Once 
again, mismatches were minimized by varying the param- 
eters of the moffat function and then performing a x^ min- 
imization. The SINFONI PSF was obtained by convolving 
the NACO PSF with the degradation function. The final 
FWHM resolution estimated in this way was 0.22" (~ 4.2 
pc) with a = 0.219" and /3 = 3.2. 

The confidence on the resulting PSF for both cases 
was evaluated by the residual of the fitting the following 
function which considers the difference between the data 
and model, weighted according the flux in each pixel, and 
normalised by the total flux in the frame: 



Residual = 



^E((IMs - Model)2 ■ IMs) 



(E(IMs)) 



1.5 



100 



where /Afs is the SINFONI continuum image as above 
and Model is our constructed image - either the intrin- 
sic galaxy proflle convolved with the SINFONI PSF, or 
the NACO image convolved with the degradation func- 
tion. The numerator is the quantity minimized during the 
fit, which was then normalized using the integrated im- 
age of the continuum. The residual of the fit in both cases 
was evaluated as 0.1%, indicating that both methods per- 
formed equally well in estimating the PSF. We chose to 
adopt the PSF obtained by the second approach because 
no deconvolution is performed. 

3.2. Distribution of emission and absorption features 

In Figure^we present a spectrum of Circinus at a spectral 
resolution of i? ~ 4200, integrated over a 0.8"square aper- 
ture centered on the nucleus. The most prominent of the 
emission lines are the coronal lines [Sivi] at 1.96 /im and 
the [Caviii] at 2.32 ^m. We also detect the [Alix] 2.04 ^m 
line, first mentioned by Maiolino et al. (1998). In addition 
to these high excitation lines, the H2 1-0 S(l) 2.12 /xm and 
Br7 2.17 /xm emission lines are clearly recognizable. The 
line fiuxes, which are presented in Tabled are comparable 
to previous measurements (Maiolino et al. ll998|l . 

The K-h&nd flux density, calculated from the contin- 
uum image shown in Figure |2] within an aperture of 
0.8" is 3.1 X 10"" Wm-2 ^.wl'^ (approximately SOmJy). 
The stellar features, traced by the CO bandheads, are 
also distinguishable in the spectrum. However, not all of 
the stellar features are diagnostically useful. For exam- 
ple, the ^^CO (3-1) band is partially filled by the [CaViii] 
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Fig. 1. Nuclear spectrum of Circinus, extracted from the SINFONI datacube in a 0.8"square aperture. Labels marked 
with an asterisk represent stellar CO absorption bandheads. The effect of dilution on these bandheads is clearly seen. 



emission line, and the rest of the bandheads with wave- 
lengths longer than 2.34 /im are affected by residual at- 
mospheric features or other emission lines. Therefore, only 
the ^^CO (2-0) band was useful to extract information 
about the fraction of the nuclear flux that is indeed stellar. 

Images of the 2.2-/xm continuum and the Br7, H2 (1- 
O)S(l), [CaVlll], and [SiVl] Hne emission, as well as the 
^^CO (2-0) bandhead flux are presented in Figure |21 The 
location of the nucleus, as defined by the centroid of the 
continuum emission, has been marked with an encircled 
cross in all maps. This source is offset by ^ 0.15" to 
the south-east of the peak of the Br7 and II2 1-0 S(l) 
line emission as well as of the stellar light. Such an offset 
is consistent with the results of Prieto et al. ( 2004 ) who 
found that the center of the i^-band radiation is located at 
^ 0.15" south-east of the brightest central emission seen 
in the J-band and HST F814W images. In contrast to 
these authors who argued that the iC-band source is fully 
obscured at shorter wavelengths and that the emission at 
these shorter wavelengths is due to nuclear light scattered 
by the compact dusty structure surrounding the Circinus 
nucleus, our results suggest a different interpretation be- 
cause the offset remains even in the K-hand but only for 
particular components of the line and continuum emission. 
We argue that, at least for near infrared emission, the off- 
set arises because the star formation is centered ~0.15" 
to the north west of the non-stellar continuum peak. The 
issue of extinction is addressed further in Section [3 



4. Nuclear dust emission 

The observed equivalent width of the ^^CO (2-0) band- 
head was used to obtain the fraction of the nuclear flux 
from the stars in this region, by assuming that the in- 
trinsic (i.e. purely stellar) equivalent width in the nuclear 
region must be consistent with values predicted by the 
models created with the stellar population synthesis code 
STARS ( Sternberg 1998| |Thornley et al. 2000[ Davies et 
al. .2003) . We used the measured flux of the ^^CO (2-0) 
band to estimate the amount of non-stellar continuum us- 
ing the equation: 



l-D 






intr 



where D is the dilution factor, that is the ratio of stellar 
to total continuum, and EWohs and E'Wintr are the ob- 
served and intrinsic stellar equivalent widths respectively. 
The intrinsic equivalent width is nearly independent of 
star formation history for any ensemble of stars as pre- 
dicted by the STARS models, having an almost constant 
value of '^ 12 A assuming solar metallicity (Davies et al. 
UOO6 ) . The observed stellar equivalent width in the central 
0.8" was found to be ~1.8 A, implying that ~85% of the 
total nuclear K-hand luminosity in this region originates 
from the non-stellar component, a fraction that increases 
on smaller scales. In Figure the K-hand continuum az- 
imuthally averaged radial surface brightness profile is plot- 
ted along with those of the PSF and stellar and non-stellar 
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Fig. 2. Intensity images extracted from the SINFONI data cube in the central arcsec of Circinus. In each case, the 
colour scale extends from 0-100% of the peak flux, and contours are spaced equally between 20% and 90% of the peak 
flux. An encircled cross indicates in each case the peak of the continuum emission. The maps show, starting from 
the left corner. Top left: 2.2/im continuum. Top center: stellar component (flux contained in the stellar absorption 
bandhead ^^CO(2-0)), Top right: Br7, Bottom left: H2 l-OS(l), Bottom center: [Caviii], and Bottom right: [Sivi]. 



components. The stellar component represents the radial 
profile of the stellar absorption bandhead ^^CO(2-0) im- 
age in Fig. [5] All profiles are normalized to their respective 
peak values; and the stellar profile is centered 0.15" north 
west of the non-stellar profile. The predominance of the 
non-stellar component on these scales is vividly demon- 
strated by the similarity of the non-stellar and total K- 
band profiles, as well as by the red spectral slope in Fig.^ 

Such nuclear ii'-band emission is interpreted as ther- 
mal radiation from hot dust which is heated by the intense 
UV and X-ray radiation emitted by the Seyfert nucleus. It 
is often assumed that the dust is close to the sublimation 
limit (~ 1600 K) since it is the hottest emission which will 
have the greatest impact in the near infrared. However, 
the slope of the nuclear spectrum in the X-band indi- 
cates a temperature (unreddened) closer to 740 K, similar 
to that found for NGC 1068 by Thatte et al. (.1997), but 
significantly more than the 300 K derived by Prieto et al. 
H2004|l from modelling the broad band colours. Since the 
spectral range of the iC-band is smaller than the 1-10/im 
interval used by Prieto et al. I)20()4|l . the spectral slopes 



in these ranges are different, but nevertheless compatible. 
The estimated temperature of ^700 K from the K-hand 



continuum characterizes only the hottest dust, while 300 K 
dust would dominate the mid- infrared SED. 

As can be seen from Figure O the non-stellar X-band 
source is marginally resolved. Indeed, a quadrature cor- 
rection of its FWHM (- 0.3") with that of the PSF yields 
an intrinsic size of ~ 0.2" (3.9 pc). The factor of 2 dis- 
crepancy with the size reported by Prieto et al. H2004|l 
arises primarily from the spatial resolution estimated in 
the NACO data. Compared to their 0.16" our estimate 
is 0.14", which would give a source size of ^ 0.15". This 
emphasises the difficulty of estimating the intrinsic size of 
only marginally resolved sources. Within the uncertainty 
of about 0.05", the sizes are consistent with each other and 
with those predicted by the torus models (Pier & Krolik 
\IM^ [BMI Nenkova et al. IM)^ in which the inner edge 
of the dust torus lies ~1 pc away from the AGN. 

5. Nuclear star formation 

Images of the Br7 and II2 1-0 S(l) line emission, which 
are often associated with star formation, are presented in 
Figure |21 As the AGN is highly obscured (see Oliva et 
al. . 1998. and Matt et al. ^19961 IT^^ . no broad line re- 
gion is visible. Several lines of evidence suggest that the 
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Fig. 3. Radial profile of the nuclear ii'-band surface 
brightness obtained from an integrated image of the 
SINFONI datacube (crosses connected by a solid hue). 
The SINFONI i^-band stellar component (same as stel- 
lar continuum) is represented by asterisks connected by 
a dotted line. The triple-dotted-dashed line represents a 
model fitted to the stellar component data, which is the 
exponential profile convolved with the PSF. Open trian- 
gles connected by a dotted-dashed line indicate the frac- 
tion of the JT-band radiation which does not come from 
stars. The open diamonds connected by a short-dashed 
line indicate the PSF profile. 



Br7 emission is associated with star formation activity 
surrounding the Seyfert nucleus rather than the narrow 
line region and ionisation cone. These are the similarity of 
the morphologies of the Br7, ^^CO (2-0), and H2 1-0 S(l) 
emission: the "'^^00(2-0) is only slightly more extended 
than the two lines, and all three are offset to the north 
west. Additional support is provided by the consistency 
and uniformity of the Br7 and H2l-0S(1) velocity fields 
and dispersion maps (see Section El and Figure [7|). In par- 
ticular, the velocity field of the Br7 has a gradient match- 
ing that of the galaxy's major axis rather than showing 
signs of outflow along the minor axis. We note in fact that 
the H2 1-0 S(l) appears to be slightly more extended than 
the Br7, particularly at weaker levels westwards in the di- 
rection of the ionization cone (Maiolino et al. 11998(1 : and it 
is only in this region that there are significant differences 
between their velocity fields, with the 1-0 S(l) showing 
signs of outflow. Figure 0] shows the azimuthally averaged 
surface brightness radial profiles of the H2 1-0 S(l) and 
narrow Br7 emission lines and the "'^^CO (2-0) stellar ab- 
sorption bandhead. 

Both an exponential profile and r^'^ de Vaucouleurs 
profile provide good matches to each of these. We prefer 
to parameterise the profile with the former, because of the 
clear evidence of a circumnuclear disk from the existence 
of certain features (namely a bar, Maiolino et al.'2000; and 
a ring, Marconi et al. ll994bjl . At these scales, less than 20 
pc, the Br7, and H2 1-0 S(l) data are characterised by a 



Fig. 4. Radial profile of the surface brightness of the H2 
1-0 S(l) (crosses connected by a triple-dotted-dashed hue) 
emission line. The long-dashed line represents a model fit- 
ted to the linemap, which is the exponential profile con- 
volved with the PSF. The narrow Br7 radial profile is 
represented by open triangles connected by a short-dashed 
line and its model by the dotted-dashed line. Open squares 
connected by a solid line indicate the nuclear stellar pro- 
file, whereas its model is represented by the dotted line. 



disk-scale length r^ = 41q'2 pc. Due to the symmetry in 
the morphology of the Br7, and H2 1-0 S(l) profiles, the 
profiles were not corrected for inclination. This would per- 
haps be a surprising result for an inclined disk, but indeed 
the kinematics and simulations suggest that the thickness 
of the disk plays an important part in the observed mor- 
phology: a disk inclined at 65° with a thickness of only 
4pc FWHM - equivalent to an exponential scale height 
of 1.7 pc - is consistent with the symmetrical morphology. 
We return to this point in Section where we discuss the 
kinematics in more detail. 

The complex morphology of the stellar component was 
analyzed in more detail by several simulations. Since the 
observed stellar luminosity is expected to be produced by 
several clusters, simulations were carried out in order to 
study the way the clusters superpose to create the ob- 
served complex morphology. In the simulations each clus- 
ter was defined by two parameters: its position inside the 
field, and its mass. The positions were assigned randomly 
following the radial distribution of the stellar component 
in Fig. 01 so that the probability of finding a cluster in 
any position depends on an exponential profile slightly 
more extended than the line emission. For the mass dis- 
tribution, systems of young clusters, including super star 
clusters which are preferentially found at the very heart of 
starbursts, appear to be well represented by a power law 
mass and luminosity function with a = —2. The range 
of masses was from 1 x lO'' to 3 x 10^ Mq which are 
typical values for star clusters (Meurer et al. 1995). Each 
cluster was characterized by a two dimensional Gaussian 
function with FWHM — 0.2", our spatial resolution. The 
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Fig. 5. Intensity images of the star clusters simulations in the central arcsec of Circinus. In each case, the colour scale 
extends from 0-100% of the peak flux, and contours are spaced equally between 10% and 90% of the peak flux. An 
encircled cross indicates in each case the center of the region. The maps show, starting from the left corner, Top left: 
1000 star clusters with constant mass 1 x 10'^ Mq, Top right: 100 star clusters with constant mass 1 x 10"* Mq, Bottom 
left: 10 star clusters with constant mass 1 x 10^ Mq, and Bottom right: 300 star clusters with masses varying from 
1 X 10^ Mq to 5 X 10* Mq. 



clusters were generated randomly over a circular field of 
0.5" radius. According to the mass distribution, a mass 
was assigned to each of them until the observed stellar 
mass in the region was reached, namely 1 x 10® Mq (see 
below). Each cluster had the same mass-to-light ratio. The 
summation of all clusters over the whole field yields the 
overall morphology of the stellar component. Fig.Elshows 
the results of the simulations. It can be seen from Fig. O 
that when massive star clusters are present in the field, the 
dumpiness of the whole set increases. On the other hand, 
for a low mass distribution the field remains uniform and 
symmetric, and therefore the dumpiness decreases. 

The observed morphology in Fig. can be best repro- 
duced by a set of star clusters with a mass distribution 
ranging from 1 x 10"^ to 5 x 10* -^0: or alternatively a 
constant mass with value of 1 x 10* Mq, although other 



scenarios cannot be excluded . Our simulations indicate 
that the observed complexity and asymmetry of the mor- 
phology of the stellar component could easily arise from 
the superposition of a number of star clusters with typical 
mass 10* Mq. 

We determine the age of the star formation by means 
of two diagnostics: the intrinsic equivalent width of the 
Br7 line and the supernova rate to stellar K-hand lumi- 
nosity ratio (vsn/L^^) in the nuclear region. These pro- 
vide independent information on the star formation his- 
tory, the former representing the number of hot young 
blue supergiant stars and the latter the number of super- 
nova remnants, both normalised to the number of later 
type red giant and supergiant stars. By comparing the 
two diagnostics to models, the age of the starburst can be 
satisfactorily constrained. The EW{Bj:'y) was estimated 



F. Mueller Sanchez et al.: SINFONI Spectroscopy of the Circinus galaxy 



1000 



S, 100 i 



CQ 




1 08 1 09 

Age (years) 



lO'O 



lO-o 



^ 10-10 



10-" 



inf 


^ 1 n8 


; V ^ 1 n6 




' jV^^ 



10' 



IQB 10» 

Age (years) 



10" 




108 109 

Age (years) 



lO'O 



1000 



100 



10 



109 

108 

10^ 

108 



10' 






108 io» 

Age (years) 



lO'O 
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from the observed Br7 flux and the stellar continuum (i.e. 
corrected for dilution by subtracting the continuum con- 
tribution from hot dust), resulting in EW (Bi^) = 30 A 
in a 0.8" aperture. The supernova rate vsn was obtained 
using the empirical relation presented in Condon (1992) 
with a spectral index a = 0.1 (as measured between 5 
and 8.64 GHz), and a flux density at 8.64 GHz of ^lo = 
12 mJy, resulting in a value of i'sn — 2.21x10"* yr"^. 
The radio flux density was calculated by means of the 1" 
resolution radio image from Davies et al. H1998|l . The flux 
density at 8.64 GHz was summed over an aperture of 0.8", 
to be consistent with our fluxes, giving 14 mJy. Then the 
thermal contribution was calculated as given in Condon 
H1992|l for a standard electron temperature of 1x10* K. 
The resulting 2 mJy was subtracted from the total flux 
density at 8.64 GHz. Dividing vsn found in this way by 
the stellar isT-band luminosity L^ of 1.5 x 10^ Lq found 
after multiplying the observed X-band luminosity Lk by 



"0 



the dilution factor D, yielded a ratio of 1.5x10 ^° 
yr"*. 

We have used the evolutionary synthesis code STARS 
( [Sternberg 1998| [Thornley et al. 200'0| Davies et al. . 2003 ,1 
to determine the star formation history assuming an 
exponentially decaying star formation rate of the form 
exp{—t/tsci), where tgci is the burst decay time scale and 
t is the age of the star cluster. We make the standard 



assumptions of solar metallicity, and a Salpeter IMF of 
slope a — —2.35 in the range 1-120 Mq. An important 
advantage of this model is that it allows us to study not 
only instantaneous and continuous scenarios, but also fi- 
nite star formation timescales. 

We obtain different ages for the starburst depending 
on the timescale of the star formation. We consider 4 
timescales: continuous star formation with tsd — 10^ yr, 
an instantaneous burst with tgci — 10^ yr, and two inter- 
mediate scenarios with tgd — 10^ and 10^ yr. First, we 
compared the calculated intrinsic Br7 equivalent width 
with that from STARS. Figure El shows the evolution of 
EW{BT'y), i^sn/L^, M^/L"^ and iboi/^k as functions of 
cluster age for each model. It can be seen from this figure 
that EW{Bi'y) alone cannot constrain the age since all 
scenarios are plausible. 

Using the ratio i'sn/L'^ together with EW(Bi'j), we 
note that within the central R < 8 pc, the constant and 
the instantaneous star formation models do not fit com- 
pletely both of the observational constraints {EW {Bi'j) = 
30 A and ^sn/^k = 1-47 x 10"!° Lq^ yr"*). 

The star formation timescales tgd = 10^ and tsd = 
10^ yr do fit these constrains from the intersections points 
at 2 x 10^ and 8 x 10^ yr respectively. While the data are 
consistent with these two scenarios, the shorter timescale 
would require a degree of fine-tuning, particularly since 
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the supernovae responsible for the radio conntinuum are 
only just beginning to appear and EW{Bij) is decaying 
relatively quickly. Thus there is only a small range of pos- 
sible ages for this timescale. At an age of 80 Myr, the con- 
straints give rise to a wider range of allowed ages, sug- 
gesting this is the more likely scenario. Crucially, in either 
case the data indicate the presence of a young stellar pop- 
ulation within a few parsecs of the active nucleus which is 
now less than half as intense as when it began. 

For the preferred star formation scenario given above 
we find a current mass to iC-band luminosity ratio of 
~O.6M0/Lo, yielding a stellar mass of M* f* 1 x 10^ Mq. 
Another important parameter to analyze is the bolomet- 
ric luminosity attributable to the nuclear star formation. 
Figure El shows the evolution of the bolometric luminos- 
ity to iiT-band luminosity ratio for an ensemble of stars. 
From this Figure we obtain iboi/^K ~ 1^0. By multiply- 
ing this ratio by the observed stellar if -band luminosity 
of L^ = 1.5 X 10^ 1/0, we find that the star formation in 
this region has L^j^j = 2.3 x 10* Lq and accounts for 1.4% 
of the bolometric luminosity of the entire galaxy (iboi — 
1.7 X 10^° Lq, Maiolino et al. .19951. 

Let us now assume a screen extinction of v4v — 9 
mag, as suggested by Maiolino et al. H1998|l for the nu- 
clear region. Since Ak « Ay/lO CHowarth 11983(1 . L^ is 
~2 times that observed. This consideration has no effect 
on the EW{Br^), but it reduces the i^sn/^^k t>y a factor of 
two. This new constraint does not change the prefered star 
formation scenario, still suggesting that the star formation 
is very young. The bolometric luminosity attributable to 
the nuclear star formation in this case would be L^^^i — 
4.5 X IO^Lq. 

A perhaps more physical model than screen extinction 
is mixed extinction, where the dust and gas are uniformly 
mixed with the stars. In either case, the near infrared color 
excess is defined as 



E„^ 



H-K 



-2.5 log 



iH/h 



OH 



Ik/Iqk 



where, Ih and Ik are the observed intensities in the two 
bands, and Iqh and Iqk are the intrinsic intensities. One 
can then use either Ix/Io\ = e""^^ for the screen model 
or I\/Io\ == (1 — e^'^^)/T\ for the mixed model to de- 
rive the optical depth. With the usual wavelength scal- 
ing that TV = 5.5th = 9.7rK one can then estimate 
the optical depth in the visual. The near infrared excess 
Eh-k = 0.45-0.60 reported by Maiolino et al. (1998) 
then implies an extinction, after accounting for the galac- 
tic foreground, of Ty — 4.4-6.2 for the screen model, and 
Ty = 12-25 for the mixed model, similar to the esti- 
mate of Prieto et al. H2()()4|l . We note also that for the 
mixed model, the colour excess will saturate because if t\ 
is greater than a few, Ix/Iox ^ I/ta so that the excess 
becomes Eh-k '^ — 2. 5 log Tif/r//. The hmiting value of 
Eh-k — 0.6 is consistent with the upper end of the mea- 
sured range, suggesting that greater extinctions cannot be 
ruled out. We can, nevertheless, impose a very strong limit 
on the maximum possible extinction: the intrinsic lumi- 



nosity from this nuclear star formation cannot exceed the 
galaxy's bolometric luminosity (iboi"*^ ^ ^boi)- The lower 
limit, that L'^^^^ > Lboi/75, comes from the starburst 
model under the assumption of no extinction. The actual 
value of ^boi"'^ ^^ governed by the scaling of the starburst 
model. And this scaling depends directly on the ratio be- 
tween the intrinsic and observed fluxes in the waveband 
where it is determined. Since the ratio between the in- 
trinsic iiT-band and bolometric luminosities of the nuclear 
stellar component is fixed by the starburst model, it fol- 
lows that the intrinsic if-band stellar luminosity must be 
less than 75 times the observed K-hand stellar luminos- 
ity. Assuming a mixed model, one has Ik/Iqk ~ ^/tk, 
implying that tk ^ 75, equivalent to ry ^ 750. A more 
meaningful, although less robust, limit is found if one ex- 
cludes the circumnuclear star formation which Maiolino 
et al. I|1998|l estimated to account for ~ 10^° L0. In this 
case Lboi/iboi"^ ^ ^^ ^^"^ ^^^ maximum optical depth 
would be of order ry ~ 300. Although these limits are 
large, they are important and further discussion on the 
extinction will be presented in Section [7] 



6. Gas kinematics and mass in the nuclear region 

In Figure [7| the projected velocity maps of the II2 and 
the Br7 lines are presented, as well as the difference in 
the velocity fields of the two lines. The maps for both 
emission lines are broadly consistent, exhibiting a gradual 
increase in velocity from north to south along a position 
angle (P.A.) of - 18° for the Br7 and ~ 25° (exclud- 
ing the righthand side, see below) for the 1-0 S(l). Given 
the typical uncertainty of 5-10° in each of these, they are 
consistent. They are the same as the galaxy's major axis 
(Freeman et al. ll977|l . a strong indication of pure galaxy 
rotation, and suggesting that there are no warps in the 
galaxy down to scales at least as small as a few parsecs. 
At smaller scales (from 0.1 to ^0.4 pc), the maser emission 
traces a warped, edge-on accretion disk, which changes the 
P.A. from -29° at 0.1 pc to -56° at 0.4 pc (Greenhill et al. 
12QilJi|. Our results suggest that at larger scales (—1 pc), 
possibly at the starting point of the torus, the galaxy's 
position angle has reverted back to —25°. 

The fact that the P.A. for the H2 1-0 S(l) velocity field 
appears to be greater is due primarily to the change in 
velocity gradient in a 0.3" strip along the western edge 
of the map; excluding this region, the difference between 
the maps is statistically insignificant (see Figure |7| right 
image). The western region coincides with the extension 
of emission in the II2 1-0 S(l) flux map and seems likely 
to be associated with a molecular outflow in (or at the 
edge of) the ionisation cone, perhaps radiatively driven 
by the AGN. It is important to notice that the observed 
velocity field of the Br7 line is clearly different to the one 
expected if it was associated with the NLR of the AGN, 
which would exhibit an increase in velocity in the western 
direction along the minor axis of the galaxy, as an outflow 
in the ionisation cone. 
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Fig. 7. Velocity maps extracted from the SINFONI data cube in the central 0.8" x 0.8" of Circinus. The maps show, 
from left to right: H2 1-0 S(l), Br7 and difference in the velocity fields of the two lines For the case of H2 1-0 S(l) 
and Br7 velocity maps, the colour scale extends from [-50, 50] km s~^, and contours are spaced equally every 10 km 
s~^. For the difference in the velocity fields, the colour scale extends from [-10, 20] km s~^, and contours are spaced 
equally every 5 km s^^ The difference in the velocity maps ranges between [-10,10] km s""'^, except for the western 
edge of the map, where a difference of ^ 25 km s~^ is observed. An encircled cross indicates in each case the peak of 
the continuum emission. 



By correcting the velocity maps for inclination (i — 
65°), we find a rotation velocity of 75kms~^ at 8 pc from 
the nucleus. The velocity dispersion of the stars in the 
nuclear region is almost constant with a value of (Jobs ~ 
80kms~^. By correcting the measured dispersion for the 
instrumental resolution corresponding to CTinst ~ 32 km 
s~^, we obtain an intrinsic dispersion of a^ ~ 70kms^^, 
yielding a velocity to dispersion ratio of Vrot/fz ~ 1.1- 
This borderline ratio suggests that while there is signifi- 
cant rotation, random motions are also important - and on 
these scales perhaps dominant. Indeed, the kinetic energy 
in random motions is of order 3a^ while that in the or- 
dered rotation is V^^^ . This would imply that the distribu- 
tion of gas is rather thick - consistent with our suggestion 
in Section |31 that the symmetry of the observed isophotes 
is due to the physical height of the nuclear region. 

Further evidence that this is the case comes from the 
scale height derived if one assumes that the gas and stars 
are self gravitating, although this should only really be 
applied to a thin disk. In this case 

a'i = 27rGSzo 

where az is the velocity dispersion of the stars, E is the 
surface density of the disk, and zq is a characteristic mea- 
sure of the scale-height of the disk perpendicular to the 
plane. In the nuclear region at i? < 0.4" (8 pc), we esti- 
mate the mass of the gas and the stars residing in the disk 
to be 1.7 X 10^ Mq (see also Section^, yielding a mean 
surface density of S = 8.5 x 10'* Mq pc~^. By considering 
an intrinsic a^ — 70 km/s, we obtain a height of zq = 2pc, 
consistent with our earlier estimate based on the isophotal 
symmetry. 

We suggest that this thickness does not necessarily im- 
ply that there is a distinct spheroidal structure in the cen- 



tral few parsecs; but that we are only looking at the very 
inner region of a larger disk which extends out perhaps 
as far as circumnuclear ring. By measuring the rotation 
velocity within a truncated field of view ~ i.e. before Vmt 
reaches its asymptotic value - one would underestimate 
Vi-ot/o". It is also possible that either dynamics associated 
with the black hole, or heating from the AGN or local in- 
tense star formation, have caused the inner region of the 
disk to thicken. 

For systems where the velocity dispersion is significant, 
estimating the mass from the rotational velocity alone 
leads to an underestimate. Accounting for random mo- 
tions is not trivial, but can be approximated by M^yn = 
(Kot + 3(7^ )r/G'. Hence we obtain a dynamical mass for 
the nucleus within ±0.4" of Mdy„ = 3.5 x lO'^ Mq. 

7. Dimensioning the clumpy star forming torus in 
Circinus 

The spatial scales on which we have traced the gas and 
stars are similar to those associated with standard torus 
models (e.g. Pier & Krohk[T992, 1993J Nenkova et al. lMil 
Schartmann et al. I^UHSl Fritz et al. IMTSJI . Indeed, the 
latter two authors have modelled Circinus specifically and 
find size scale for the torus of 30 pc and 12 pc respectively. 
However, to understand better the relation between the 
gas and the stars, we need, in addition to the stellar mass 
found in Section |31 to estimate the gas mass. Fortunately, 
as we discuss below, there is a way to do this from our 
data. 

TableElhsts 17 galaxies for which there are in the liter- 
ature both measurements of the H2 1-0 S(l) line flux and 
estimates of the total gas mass from mm CO 1-0 or CO 2- 
1 luminosities. The H2 1-0 S(l) luminosities are primarily 
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Table 2. Properties of the sample of luminous and ultralmninous infrared galaxies used to estimate i(i-o)s(i)/-^^gas 
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Fig. 8. Left panel Nuclear spectrum around [Caviii] to show how the stellar continuum was substracted. Right panel 
[CaViii] profile in velocity space, fitted by two gaussians 



from Goldader et al. H1997fl . since these are taken in large 
(3" or more) apertures and therefore are more likely to 
include all, or nearly all, the H2 1-0 S(l) flux. In addition, 
these galaxies are actively star forming ~ as evidenced by 
their classification as luminous or ultraluminous galaxies. 
As the table shows, nearly all of these have similar ratios of 
the 1-0 S(l) luminosity with the infrared luminosity Lir, 
and also with the molecular mass estimated from the CO 
luminosity Lqo- The former relation has previously been 
noted by Goldader et al. (,1997,) . who suggested it could be 
explained in terms of supernova remnants if the H2 origi- 
nates in gas shocked by the expanding shells. On the other 
hand, Davies et al. (2003) argued that in general the H2 
1-0 S(l) line in such galaxies is excited by fluorescence by 
young stars of gas dense enough that the lower vibrational 



levels are thermalised. In either case, if the H2 emission is 
primarily associated with star formation, it would not be 
unreasonable to find a correlation with Lir. 

Similarly, a correlation between Ljr and Leo has been 
known for a long time (see Young & Scoville I1H91I for a 
review). Recent work (Gao & Solomon 2004) shows that 
while LiR and Leo are correlated, the logarithmic relation 
is linear only over limited luminosity ranges and there is 
a clear trend of Li^/Lco (i-e- star formation efficiency) 
with luminosity. They also showed that if one uses instead 
a tracer of dense gas such as HCN, the equivalent ratio 
Lir/Lhcn is constant; although for our application one 
then has the difficulty of estimating the gas mass from 
the Lhcn- 
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In both relations, NGC 6240 is a major exception, hav- 
ing by far the largest H2 1-0 S(l) luminosity. The reason 
is most probably that the excitation mechanism is most 
probably shocks rather than fluorescence (e.g. Sugai et al. 
I1997|l . For this reason, we have excluded NGC 6240 from 
the sample. 

Given the correlations above, it would not be unrea- 
sonable also to expect a relation between 1-0 S(l) lu- 
minosity and H2 mass as traced by the CO luminos- 
ity (or even better would be HCN luminosity). In fact 
we find that the logarithmic ratio of these two quanti- 



ties is log- 



^(l-O)S(l) 



-3.6 ± 0.32 which corresponds to 



^ 2.5 X 10~^Lq/Mq with a la uncertainty of a factor 2. 
This is a quite remarkable result for such a diverse sample 
- that includes galaxies with powerful AGN for which X- 
ray irradiation of gas may be an important contributor to 
the 1-0 S(l) luminosity. It gives us confidence that we can 
use the H2 1-0 S(l) line, without needing to know the de- 
tails of how it is excited, to make at least an approximate 
estimate of the total (cold) molecular gas mass. 

We apply this result to the nuclear region of Circinus 
because it has had recent vigorous star formation close to 
the AGN, and thus matches the sample. From the total 
H2 1-0 S(l) luminosity of 4.4 x 10"^ Lq we estimate the to- 
tal molecular mass to be 1.7 x 10^ Mq - fully consistent 
with the dynamical mass estimated in Sectional The gas 
mass surface density is then 8 x 10** Mq pc~^, yielding 
a column density of tih = 5.2 x 10^'' cm~^, which is very 
similar to that implied by X-ray observations (e.g. Matt 
et al. ll999|l . Since nn « 1.5 x 10^"'^ ry, the mean gas den- 
sity implies an optical depth through the thickness of the 
disk of TV = 3000. This is far larger than the most likely 
values estimated in Section and also inconsistent with 
the maximum possible optical depth. 

This discrepancy suggests that the canonical torus in 
the unified model of AGN is not only forming stars but is 
a clumpy medium - so that much of the star light is not 
obscured despite the presence of huge gas column densi- 
ties. 



8. Coronal lines 

Because of the high ionization potential (IP>100 eV) as- 
sociated with Coronal Line Region (CLR) emission lines, 
highly energetic processes are required. The lines can be 
excited either by a hard UV to soft X-ray continuum, very 
hot coUisionally ionized plasma, or a combination of both. 
If the excitation occurs via collisions, the gas temperature 
should be about T — 10^ K. In the case of photoionization 
via the hard AGN continuum, temperatures of only a few 
10^ - 10"^ K are needed. 

In our observations of Circinus the [Sivi] , and [Caviii] 
lines are extremely strong and [Alix] is also detected, 
as apparent in the spectrum in Figure Q Due to the 
low signal to noise ratio of the K-hand spectrum after 
2.42/zm, we exclude the [Sivii] line at 2.48 /im studied by 
Prieto et al. 1)2005(1 . From this spectrum and also that in 
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Fig. 9. Comparison of coronal line profiles in velocity 
space. 
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Fig. 10. Radial profile of the surface brightness of the 
[Caviii] broad component (asterisks connected by a dotted 
line), and [Caviii] narrow component (open triangles con- 
nected by a dashed line) , compared to the [Caviii] profile 
(crosses connected by a dotted-dashed line) and the PSF 
profile (open diamonds connected by a solid line) 



Figure IHl it can be seen that the [Caviii] line sits on top 
of the stellar ^^ C0(3-l) absorption bandhead. A template 
star (HD 179323, shifted to the redshift of Circinus) was 
used to correct the stellar features in the spectrum and 
hence reconstruct the full [Caviii] line profile. As the deep 
CO bandheads are produced in the atmospheres of red 
(super)giant stars which dominate the emission around 
2.3 /xm, HD 179323 (spectral type KOIab) was choosen for 
this purpose. The correction was achieved by convolv- 
ing the template with a Gaussian broadening function, 
and varying its parameters to minimise x^, which was 
measured across the 1^00(2-0) bandhead. The ^"^00(3- 
1) bandhead of the template was then convolved with 
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Line 


IP 


Component 




Flux 


-Abroad / -^narrow 


Velocity'' 


FWHM Line Wdith 




(eV) 




(W 


m-2xl0-^«) 






(km s^^) 


(km s"^) 


(1) 


(2) 


(3) 




(4) 




(5) 


(6) 


(7) 


[Sivi] 


167 


Broad 




26.8 






-103 


300 


[Sivi] 


167 


Narrow 




24.7 




1.1 


~0 


175 


[Mix] 


285 


Broad 




2.6 






-228 


300 


[Mix] 


285 


Narrow 




8.2 




0.3 


~0 


200 


[Caviii] 


127 


Broad 




28.0 






-87 


540 


[Caviii] 


127 


Narrow 




20.0 




1.4 


-0 


150 



" The ratio for each line is written in the row of the narrow component. 
' Velocity shift from the systemic velocity. 

Table 3. Properties of the broad and narrow components of the detected coronal lines. 



the broadening function and subtracted from the galaxy's 
spectrum. 

Figure shows that in all three coronal line profiles 
we observe asymmetric and broadened lines, indicating 
the presence of two or more components within the coro- 
nal line region. In order to quantify the profile, we have 
fitted a superposition of two Gaussians to each line, as is 
summarised in Tabled As can be seen in Figure |S1 for the 
case of [Cavili], there is a strong narrow component and 
a weaker blueshifted broad component. These also have 
different spatial extents. Figure ITUl shows azimuthally av- 
eraged radial profiles of the broad and narrow components 
of this line, as well as the PSF. It is clear that the broad 
component is resolved and spatially more extended. The 
flux peaks at about 1.5 pc from the nucleus. On the other 
hand, the narrow component is unresolved and most of 
its flux comes from a region of less than 1 pc across and 
coincident with the non-stellar K-hand continuum peak 
(rather than the peak in Br7, II2 1-0 S(l), and stellar con- 
tinuum). These characteristics suggest that they originate 
in different regions around the AGN, and perhaps are even 
excited by different mechanisms. 

In the case of the narrow component, the fact that it 
is compact and centered on the nucleus suggests that it 
originates physically close to the AGN. Since the lines are 
narrow and at systemic velocity, these clouds of ionized 
material are not out-flowing which is also an indication 
that they are excited by photoionization (Oliva at al. 119941 

In the case of the broad component, the fact that it 
is blue shifted indicates that part of the gas must arise 
in outflows around the AGN, as has been suggested by 
Rodriguez- Ardilla et al. (•2,0D.4). Previous work on the coro- 
nal line emission in Circinus (Oliva et al. I1994|l had as- 
cribed the excitation to photoionisation. This was in part 
due to the narrow width of the most prominent compo- 
nent (typically 175kms~^, from Tabled. However, the 
considerable broadening associated with the blue wing 
(FWHM>300 km s^^) opens the possibility that some 
fraction of it might in fact be excited by fast shocks. The 
picture of a shock-excited blue component is supported 
by the results of Prieto et al. (|2005|l . who suggested that 



in addition to photoionization, shocks must contribute to 
the coronal emission. They proposed shock excitation as 
an additional energy source to explain the extended coro- 
nal line emission, as traced by [Sivii], in a sample of 4 
nearby Seyfert 2 galaxies. Although it could be a single 
broad component, we consider a more likely scenario to 
involve many narrow components arising from different 
clouds moving at different velocities whose emission com- 
bines to produce the observed morphology and profile. A 
mechanism by which broad coronal lines are created in 
out-flowing cloudlets which have been eroded from the 
main clouds is an intepretation which has been proposed 
for the detailed morphological and kinematic data avail- 
able for NGC 1068 (CecU et al.EOIE)). 

It is interesting that we observe a correlation between 
the IP and the blueshift, as well as the i^broad /^narrow, of 
the forbidden lines. As the ionization potential increases, 
the blueshift also increases and the i^broadZ-Fnarrow de- 
creases. This can be understood if the acceleration of 
cloudlets away from the AGN is driven impulsively (which 
one might expect given the enormous variability of hard 
X-ray luminosity in many AGN) so that the cloudlets with 
the highest velocities will also be closest to the AGN - and 
as they travel further out, they will be decelerated by drag 
against the interstellar medium (probably small in the ion- 
isation cone) and by the increasing gravitational pull. The 
gravitational pull is signiflcant since it depends not only on 
the black hole mass, as assumed by Maiolino et al. H2()()()|l . 
but also on the gas and stars in the nuclear region which 
we have shown provide an order of magnitude more mass 
than the black hole alone on scales of 10 pc. Additionally, 
those lines which required the hardest ionising continuum 
will arise preferentially closer to the AGN. Thus although 
the strength of the narrow component which originates 
close to the AGN will depend on IP since there are fewer 
photons energetic enough to ionise species with higher IP; 
it will be even harder to create these lines in outflowing 
cloudlets, since they are already further from the AGN: 
they will arise only in the fastest clouds closest to the 
AGN. Lines with less extreme IP will be generated also in 
slower clouds further from the AGN and hence the mean 
blueshifted velocity will be less. Although we do not find 
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a clear correlation between the IP and the width of the fit 
to the broad component, that of the [Caviii] line exhibits 
the highest dispersion for the lowest IP, providing further 
support for this hypothesis that lines with the lower IP 
can arise in regions which are further from the AGN and 
hence moving outward slower. 

9. Conclusions 

We have presented near-infrared adaptive optics integral 
field spectroscopy of the nuclear region of the Circinus 
galaxy with an angular resolution of 0.2" and spectral 
resolution R ^ 4200, which we use to probe the gas and 
stellar morphologies and kinematics on scales of a few par- 
sees. 

In the central 0.8", the non-stellar continuum domi- 
nates the X-band, contributing 85% of the total flux den- 
sity. Offset by -0.15" from its peak are the H2 1-0 S(l) 
and Br7 lines and stellar continuum. The similarity of 
their morphologies and kinematics lead us to conclude 
that they all originate predominantly with recent vigor- 
ous star formation. Analysis of these diagnostics together 
with the radio continuum suggest that within 8 pc of the 
AGN there is a starburst which is exponentially decaying 
on a timescale of 100 Myr and began only 80 Myr ago, and 
currently has a bolometric luminosity of —1.4% that of the 
entire galaxy. 

We have found that it is possible to make a reasonable 
estimate of the total gas mass from the H2 1-0 S(l) line 
luminosity, finding that it is an order of magnitude more 
than the young stars, and also consistent with the column 
densities to the AGN implied by X-ray observations. Since 
the scales on which the gas and stars exist arc similar, we 
suggest that the torus is forming stars; and because of the 
limits on extinction to the stars, that the torus must be a 
clumpy medium rather than uniform. 

The coronal lines comprise two components: a promi- 
nent narrow (FWHM ~ 175kms~^) part which is at sys- 
temic velocity, spatially unresolved (at our 4pc FWHM 
resolution), and centered on the non-stellar continuum; 
and a broad (FWHM > 300 km s^^) blue shifted part 
which is spatially extended. We argue that the narrow 
part arises in clouds physically close to the AGN; and 
that the blue wing originates, as appears to be the case in 
some other Seyfert galaxies, from cloudlets that have been 
eroded form the main clouds and are accelerated outward. 
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